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Abstract 

In this note, we develop a theory of Euler-Poincare reduction for discrete 
Lagrangian field theories. We introduce the concept of Euler-Poincare equations 
for discrete field theories, as well as a natural extension of the Moser-Veselov 
scheme, and show that both are equivalent. The resulting discrete field equations 
are interpreted in terms of discrete differential geometry. An application to the 
theory of discrete harmonic mappings is also briefly discussed. 



MSG classification (2000): 58H05; 65P10; 70S10 



1 



1 Introduction 



Over the past decades, the interest in discrete mechanics has been steadily increasing. 
Partly this is due to the development of geometric integrators, which provide a relatively 
simple means for the long-term numerical integration of mechanical systems. On the 
other hand, many concepts (Poisson tensors, reduction, etc.) from the continuum theory 
also have a natural counterpart in the discrete realm, which makes discrete mechanics a 
subject worthy of interest in its own right. We refer here in particular to the pioneering 
efforts of Moser and Veselov [17], who studied a class of discrete integrable systems, 
and Weinstein [21], who was the first to recognize the fact that many of these discrete 
systems can be understood from the point of view of Lie groupoids. 

Inspired by these powerful and elegant methods, a number of people set out to develop 
a similar geometric approach to classical field theories. Bridges [3] introduced a concept 
of "multisymplecticity" for Hamiltonian partial differential equations and later Bridges 
& Reich [4] studied numerical integrators that conserve a discretized version of multi- 
symplecticity. Independently, Marsden, Patrick, and Shkoller [TB] extended the work 
of Veselov in order to deal with Lagrangian field theories. In a previous paper [20], we 
built upon their work by extending the Lie-groupoid methods from Weinstein's paper 
and related work by Marrero, Martin, and Martinez |14j to the case of Lagrangian field 
theories, allowing us to treat a larger class of field theories, as well as shedding new 
light on some of the constructions in [TB] . 

In this note, we now focus on a special class of discrete Lagrangian field theories and 
their behaviour in the presence of symmetry: we introduce a reduction procedure, 
similar to the one proposed by Marsden, Pekarsky, and Shkoller [TB] and Bobenko and 
Suris [H [2] for discrete mechanical systems, and show that the resulting discrete field 
equations have a simple and natural interpretation in the context of discrete geometry. 

We begin by introducing a fixed mesh (a certain collection of vertices, edges, and faces) 
in the space of independent variables. Initially, we are interested in discrete fields 
that associate to each vertex an element of a given Lie group G. If the Lagrangian 
of such a model is G- invariant, we show in ^ that these fields can be reduced to a 
new class of fields that associate a group element to each edge of the given mesh. This 
reduction procedure is a field-theoretic version of the discrete Euler-Poincare equations 
from [15\. We then show that the converse procedure (of reconstruction) is possible 
only if a certain obstruction vanishes. This obstruction has a natural interpretation as 
the curvature of a discrete G-connection. This is similar to the case of continuum field 
theories, as studied in [3 [6] . 

In §H we take an alternative route to the Euler-Poincare equations: inspired by a 
similar treatment in [5], we prove a discrete version of Noether's theorem and use the 
G-invariance of the Lagrangian to derive the field equations. Finally, in ^ we consider 
the Lagrangian of harmonic mappings into a Lie group: we propose an extension of the 
well-known Moser- Veselov algorithm, demonstrate its equivalence to the Euler-Poincare 
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equations, and establish a particularly clear form of the field equations (involving the 
concepts of discrete geometry introduced in 

Throughout the paper, a number of simplifying assumptions will be made. In ^ we 
therefore conclude this paper by giving a brief overview of how these assumptions may 
be circumvented, and we also indicate how the discrete field theories in this paper fit 
into the framework of Lie groupoid field theories developed in [20]. In particular, we 
show that the Euler-Poincare reduction procedure in ^is just a particular instance of 
a general reduction theorem in the category of Lie groupoid field theories. 

2 Discrete differential geometry 

2.1 Discretizing the base space 

In this paper, we consider field theories with two independent variables. Such field 
theories can be modeled as sections of a fibre bundle with base space (see [T6]). 
In order to discretize these fields, we introduce the set of vertices V in as the 
set of points in with integer coordinates {V = Z x Z). As in finite-difference 
approximations, the idea is that continuous fields are approximated by specifying their 
values on the elements of V. 

For future reference, we also introduce the set of edges E, whose elements are ordered 
pairs of the following form: for all i, j G Z x Z, 

+ e E, and + 1)) e E. 

Furthermore, we also demand that if {xo,xi) G E, then also (a;i,xo) G E, and we 
write (xi,xo) = {xo,xi)~^. In other words, E consists of "horizontal" and "vertical" 
line segments of unit length. To avoid cumbersome notation, we will note an arbitrary 
element of E as e, or if we have to refer to its begin and end vertex, as {xq, Xi). 

Finally, we introduce the set of faces F as the set of quadruples (cq, Ci, t2, es) of closed 
paths in E, i.e. where the end vertex of is the begin vertex of ej+i for z = 0, . . . , 3, 
and, in addition, we demand that ^ e,^\ for z = 0, . . . , 3. (The indices should be 
interpreted as being "modulo 4": C4 is just eo.) 

The elements of F can thus be pictured as little squares in M^. A generic element 
of F will be noted as f, or, since f is fully determined by its bounding vertices, as 
{xo,Xi,X2,Xs). Note that there are two distinct classes of elements in F: those where 
the vertices are denoted in clockwise and anticlockwise fashion, respectively. 

2.2 Discrete differential forms 

This section is devoted to a review of some elementary concepts of algebraic topology, 
which will be needed in our discussion of discrete harmonic maps in section [51 Most of 
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this section is modeled on [8], as well as on the text book [9]. 

Consider the mesh (V, E) in introduced in the previous section. The collection of 
sets {V, E, F} together with its various incidence relations determines a cell complex 
and this leads us naturally to the concepts of homology and cohomology. It is therefore 
not unreasonable to expect that some of these concepts will enter our study of discrete 
field theories later on. 

An n-chain on the given mesh, where n = 0, 1, 2, is a formal linear combination (with 
coefficients in R) of "n-dimensional elements". More precisely, the vector space of 
0-chains is given by 

Co = {aixi H h araXm : ai, . . . , a„i e M, xi, . . . , e V], (1) 

where it should be stressed that the elements of Co are formal linear combinations of 
elements of V . Similarly, the vector space Ci of 1-chains is generated by elements of 
C2 by elements of F, and so on. Furthermore, for the space of 1-chains, we adopt 
the identification that — (xo,xi) = (xi,xo). 

It is customary to define discrete n-forms as ra-dimensional cochains, i.e. elements of 
the dual vector space C* of C„. From this definition, it follows immediately that a 
discrete zero-form induces a function (p : V —>■ M.. Conversely, such a function gives rise 
to a zero-form through linear extension. 

Similarly, discrete one-forms can be identified with functions (p defined on the set of 
edges. It should be borne in mind that these functions are not necessarily defined on 
the whole of x V^, but only on the subset E of edges of the mesh. We also note that 
V?(xo,xi) = -ip{xi,xo) for any edge (xo,xi). 

We continue by defining discrete two-forms as functions ■?/'(a;o, Xi, X2, ^3) on the set of 
faces. Again, these functions can be extended unambiguously by linearity to a proper 
cochain. To summarize, we have the following definition: 

Definition 2.1. For n = 0, 1, 2, a discrete n-form is a linear map / : C„ — > M. For 

n > 2, all discrete n-forms are zero. The set of all discrete n-forms is denoted by C*. 

For the sake of self-containedness, we recall the explicit form of the coboundary operator 
d : C* — *• C*^i. For a zero-form 0, d0(xo, Xi) = (p{xi) — 4>{xo). For a one-form ip, 

d(p{f) = Lp{xo, Xi) + (p{Xi, X2) + 'P>{X2, X3) + Lp{x3, Xq) for all f = (Xo, Xi, X2, X3). (2) 

The coboundary of a discrete two-form is defined to be zero. We will sometimes refer 
to d as the "discrete differential". Note that d is the dual of the boundary operator 
d : C„+i Cn, in the sense that {df,v) = {f,dv) for all f E C* and v G C„, and 
where (■, ■) is the natural pairing between Cn and C*. This is a discrete analogue of the 
Stokes theorem. For more information, see [9l p. 186]. 
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Figure 1: Square mesh (black) and its dual (light). 



2.3 The discrete Hodge star 

In Riemannian geometry, the Hodge star on an oriented manifold M maps n-forms 
on M to (m — n)-forms on M, where m = dimM. As could be expected from the 
continuous theory, the discrete Hodge star to be introduced below, maps discrete 
n-forms into (2 — n)-forms. However, there is an additional complication in the discrete 
case: the forms -kf are not defined on the mesh itself, but rather on a dual mesh, which 
we now define. 

The dual mesh {V*,E*) is constructed as follows. For every face in F, there is a 
vertex in V* (for which one usually takes the circumcentric dual; see |8j). There is an 
edge in E* between two vertices qo,qi G V* if and only if the faces in F corresponding 
to go and qi have precisely one edge in common. This determines the sets V* and E*; 
the set F* consists of the faces of this dual graph. It is easy to see that to each face in 
F*, there corresponds a vertex in V. See figure [1] for an illustration. 

Implicit in this definition is the existence of a duality operator * between n-chains on 
the mesh and (2 — n)-chains on the dual mesh. This duality is well defined, but only 
up to orientation; we now use the standard orientation of to settle this point. The 
definition used here agrees with the algorithm for the orientation of dual cells proposed 
in [ini remark 2.5.1]. 

We first define * on the elements of V, E, and F. By linear extension, it will then 
be determined on the whole of Cn- The dual vertex *f of a face f whose vertices are 
written down in anticlockwise fashion is the element r of V* introduced above. If, on 
the other hand, the orientation of f is the opposite, then *f equals — r. We define the 
dual *x of a vertex x to be the corresponding face in F*, with the natural orientation: 
*x = (ri, . . . ,r4). Finally, the definition of * on is slightly more intricate; here, we 
follow [To]. Let {xo,Xi) be an edge in E and let {ro,ri} be the corresponding dual 
edge, considered as an unordered set. The line segments [a;o,a;i] and [ro,ri] determine 
a basis of M?: if this basis is positively oriented, then *{xo,Xi) = (ro,ri), otherwise, 
*{xo,xi) = (ri,ro). In the case of the square mesh of figure [H, the action of * on 
corresponds to an anticlockwise rotation over 7r/2. 

On the dual mesh, one can again introduce discrete forms. We will denote the vector 
space of discrete n- forms on the dual mesh by D*. As the dual mesh is again a square 
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mesh, there is a natural way to extend * to an operator from {V*, E*) to (V, E). It is 
then easy to check that * * v = (^—l^^(^~"'^v for any v E Cn- 

Definition 2.2. The discrete Hodge star : C* ^ ^■^ defined by 

{•ka){*v) = a{v) for all v G C„. 

The definition given here is (up to a constant) a special case of the one proposed in [8]. 
Note that *i.a = 

With the discrete Hodge star and the coboundary operator of the previous paragraph, 
we now arrive at the definition of the discrete codifferential. 

Definition 2.3. Let a be a discrete n-form. Then the discrete codifferential 6a is 

the discrete {n — l)-form 6a defined as 6a = -kd-k a. 

It is useful to write out this definition explicitly for one-forms and two-forms. If is a 
discrete one-form, then 6ip is given by 

{6lp){x) = ip{Xi,x) + (f{x2, X) + ^^(Xs, X) + V3(X4, x) , 

where Xi,X2,Xs, X4 are the end points of the edges that emanate from x. In other words, 
6ip assigns to each vertex x the sum of contributions from on the edges that have x as 
a vertex. Secondly, for a discrete two-form we note that 6iIj is the discrete one-form 
given by 6ip{xo,xi) = ip{fo) — 'ipifi), where fo and fi are the faces that have {xo,xi) as 
a common edge, and where fo is the face where the orientation of the boundary edges 
agrees with the ordering of {xo,xi), whereas fi is the face with the opposite ordering. 

Remark 2.4. Our definition of discrete codifferential agrees with the one from Desbrun 
et al. ^ up to sign. In their paper, the discrete codifferential is defined as 6a = 
(— l)*^""*"^ -k d-ka, where k is the dimension of the ambient space. Note that here, as 
k = 2, (— = —1 regardless of n. o 

2.4 Discrete connections 

In the preceding sections, we introduced discrete one-forms gnments of a real 

number to each edge t E E. This theory can be extended in a straightforward way to 
discrete forms taking values in an arbitrary Abelian Lie group G, the only significant 
difference being that we have to redefine the spaces of n-chains C„ as consisting of 
formal linear combinations with coefficients in Z. 

For instance, if <y9 : Ci ^ G is a discrete one-form, then dy? is determined by its action 
on the set -F by ([2]) and can be extended by linearity to yield a map from C2 to G, where 
it should be borne in mind that the elements of C2 are still formal linear combinations 
of elements in F, but now with coefficients in Z. 
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The theory of discrete forms with values in an Abehan Lie group will be used in section[5], 
but in the general case, we will be confronted with mappings from ii^ to a non-Abelian 
Lie group G. Such maps can no longer be interpreted as discrete one-forms. Luckily, 
it turns out that these maps have a natural interpretation as discrete G -connections, 
which we now define. 

Definition 2.5. A discrete G-connection is a map uj : E ^ G, such that, for 
all edges t E E, ^^(e"^) = uj{t)^^ . The curvature of such a connection is the map 
Q : F ^ G defined as f2(f) = Lj{ti) ■ ■ ■ u{ti), where Ci, . . . , 64 are the boundary edges of 
the face f. A discrete G-connection is said to be flat if = e for all f € -F. 

Note that in the case of a non-flat connection, f2(f) depends not only on f, but also 
on the exact representation of f as a set of edges ei, . . . , 64 (any cyclic permutation of 
this set represents the same face). However, this indeterminacy does not occur for fiat 
connections, the only case that we will consider later on. 

The theory of discrete G-connections closely mimics the usual theory of connections. As 
an example, we mention the following proposition, from which a number of interesting 
properties may be deduced. 

Proposition 2.6. Consider a discrete G-connection u : E G. If uj is flat, then there 
exists a unique mapping cf) -.V xV ^ G such that (j)\E = uj- 

This follows immediately from [20l prop. 7], or can be proved directly as follows. 

We define a path in to be a sequence Ci, e2, . . . , Cm of edges, such that the end vertex 
of Ci is the begin vertex of e^+i (for i = l,...,m — 1). Let u : E ^ G he any discrete 
G-connection. In particular, uj need not be fiat. Then, u induces a map u from the set 
of paths to G as follows: 

(£;(ci, e2, . . . , tm) = uj{ti)uj{t2) ■ ■ ■u}{tm). 

We call u the "discrete holonomy mapping". Note that if u; is a fiat connection, one 
can easily prove that Co maps closed paths to the unit in G. Furthermore, in the case 
of a fiat connection, it can be easily established that simplicially homotopic paths have 
the same image under uj. (This is the discrete counterpart of a well-known theorem of 
continuous connections (see [UJ p. 93]): if a; is a fiat connection, then any two closed 
homotopic loops have the same holonomy.) 

Let us now return to the discrete theory. In the case of a fiat connection on a simply- 
connected cell complex, we have seen that uj does not depend on the choice of path 
between two fixed points, and induces therefore a map uj : V x V ^ G. This map 
satisfies the requirements in proposition 12. 6[ 

Remark 2.7. The concept of discrete G-connections used here is common in lattice 
gauge theories (see [22] )• Novikov [T8] studied a similar concept. 
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Closely related to this definition is the concept of discrete connection on a discrete 
principal fibre bundle (see [12] for definitions and [TH §5.4] for an application to discrete 
reduction). The latter could be used to extend the results in this paper to the case of 
Lagrange-Poincare reduction, where a G-invariant Lagrangian on a pair groupoid QxQ 
is given. If Q coincides with the symmetry group G, then Lagrange-Poincare reduction 
is just Euler-Poincare reduction, the case considered here. <> 

3 Discrete Euler-Poincare reduction 

In this section, we begin our study of discrete Lagrangian field theories. Initially, we 
define discrete fields as follows: 

Definition 3.1. A discrete field is a map (p : V ^ G. 

If the discrete Lagrangian, to be defined below, is G-invariant, we shall see that these 
fields induce a new class of discrete fields, that associate a group element to each edge. 
From the last section, we know that such maps have a natural interpretation as discrete 
G-connections. 

In theorem 13.41 it is shown how the field equations for the unreduced fields are 
equivalent to a set of equations, called discrete Euler-Poincare equations, for the reduced 
fields (f. Both sets of equations arise by extremizing a certain action functional. In 
theorem 13. 6[ we deal with the reconstruction problem. Starting from a reduced field 
ip : E ^ G, it is shown that ip gives rise to a solution (p : V ^ G of the original field 
equations if and only if the curvature of if vanishes. This treatment was inspired by the 
work of Castrillon and Ratiu [6] , who developed Lagrangian reduction for field theories 
in the continuous case. 

3.1 The discrete Lagrangian 

Definition 3.2. A discrete Lagrangian is a function L : G x G x G ^ 

A word of explanation is in order here. In the continuous Lagrangian is a 

function C : J^ir — > M, where J^vr is the first jet bundle of the bundle ti : x G M^. 
It can be shown that, in this case, J^vr is isomorphic to x [TG © TG) (see p!9l 
lemma 4.1.20]). Under the (modest) assumption that Lagrangian C does not depend 
on the coordinates of the base space, we therefore conclude that C is just a function on 
TG © TG. 

The idea of defining a discrete Lagrangian as a function on G x G x G then follows 
from the idea of Moser and Veselov of approximating TG hj G x G: by applying this 
Veselov-type discretization twice (once in the "horizontal" and once in the "vertical" 
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Figure 2: Discretization of TG ® TG by G x G x G 



direction of the square mesh), it follows that TG © TG has G x G x G as a natural 
discrete counterpart; see figure [2l 

The Lie group G has a natural diagonal action by left translations on the Cartesian 
product G X G: g ■ ((71, (72) = {991,992)- We now wish to study the situation where L is 
invariant under this action. As the quotient {G x G)/G is naturally isomorphic to the 
Lie group G itself, a discrete field (f) : V ^ G induces a "reduced field" ip : E G as 
follows: 

(p{t) = (l){xo)~'^(j){xi), for all e = {xq, xi) G E. 
This leads us to the following definition: 

Definition 3.3. A reduced discrete field is a map (f : E ^ G such that f{t^^) = 
ip{t)^^. Equivalently, it is a discrete G-connection. 

The reduced field if induced by a discrete field (p -.V ^ G is fiat when considered as a 
discrete G-connection. Hence, by proposition 12.61 (or simply by construction), it may 
be extended to a map ip -.V xV ^ G. However, not all reduced fields are fiat. 

The Lagrangian L gives rise to a reduced Lagrangian / : G x G ^ M defined by 



Let $ : G X G — s> G be the map defined as ^{9, 9') = 9 ^9', so that the reduced field ip 
associated to a discrete field is given by 



l{9i ^92, 9i ^93) = L{gi,g2,93)- 



$ o 



(0 X 0). 



(3) 



Furthermore, we introduce the following mapping: 



^ : {gi,g2,g3) eGxGxG^ 



{9i 92, 9i gs) eGx G. 



(4) 



It is obvious that the Lagrangian L is related to the reduced Lagrangian by L 
For future reference, we also introduce a map ^ : T(G^) —>■ {G x g)^ as 



'^{vi,V2,V3) = {gi,g2,g3;TLg-i{vi),TLg-i{v2),TLg-i{v3)), 
where Vi G Tg.G. 



(5) 
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An unreduced field (p : V ^ G can be interpreted as an assignment of an element (pij 
in G to eacli vertex (i, j). Similarly, a reduced field ip : E ^ G can be described as an 
assignment of a group element Uij to each "vertical" edge ((z, j + 1), and of a 

group element Vij to each "horizontal" edge ((i, j), {i + 1, j)), where 

"ij = (Pij^Pi+iJ and = (6) 
We will use these notations for the remainder of the paper. 



3.2 The reduction problem 

Given a discrete Lagrangian L, the discrete action sum S is given by 

^i^) = X] L{(pi,j,(pi+i,j,(pi,j+i), 

where is a map from V to G, and where [/ is a finite subset of the set of vertices V . 
In [1^, it was shown that (j) is an extremum of this action if and only (f) satisfies the 
following set of discrete Euler-Lagrange equations: for all (z, j) G V, 

dL dL dL 

Similarly, we may define the reduced action sum s as 

A reduced field : E ^ G is an extremum of s if and only if it satisfies the discrete 
Euler-Poincare equations, to be derived below. The central aspects of discrete Euler- 
Poincare reduction are summarized in the following theorem. This theorem, as well as 
its proof, are very similar to the discrete reduction process in mechanics (see |1 5J ) . 

Theorem 3.4 (Reduction). Let L he a G-invariant Lagrangian on G x G x G and 
consider the reduced Lagrangian I on G x G. Consider a discrete field (p : V ^ G and 
let (f : V X V ^ G be the induced reduced field. Then the following are equivalent: 

(a) (j) is a solution of the discrete Euler-Lagrange equations for L; 

(b) (j) is an extremum of the action sum S for arbitrary variations; 

(c) the reduced field (f is a solution of the discrete Euler-Poincare equations: 



(/?;//(-, t;„))^-(L;_^//(-,t;,_i,,) 
(i?;//(n,„-))^-(L;,^.,d/(u 



+ 

= 0; 
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(d) the reduced field (p is an extremum of the reduced action sum s for variations of 
the form 

5ui^j = TLu,^^ {Oi+.j) - TRu,^^ {Oi,,) G T„,_^.G (9) 

and 

6v,,, = TL,^^^ (^,,,+1) - TR,^ ^ (^,,,) G r,,_^.G, (10) 
where dij = TL^-i{6(j)ij) G g. 

Proof. The equivalence of ^ and (jb]) follows from a standard argument in discrete 
Lagrangian field theories, and was shown in [20] • 

In order to prove the equivalence of (jb]) and ([d]), we note that L = lo^ (see (jl])), from 
which we conclude that if = $ o 0, then 5(0) = s{(p). Now, consider the components 
{ui^j} and {vjj} of the reduced field, as in ([H]). It is easy to check that an arbitrary- 
variation e I— >■ (^ij(e) of (f) induces corresponding variations and 5vij oiuij and Vij, 
given by fl9D and (flUj) . 

Finally, we show the equivalence of ([d]) with the Euler-Poincare equations ([H]). From 
the reduced action sum s{ip) we obtain that 



d 



d 



€=0 ^ de 



e=0 



(d/(-, T;ij) ■ bui^j + d/(Mij, ■) ■ bvi^j) . 



Substitution of ^ and (ITU]) into this expression then yields (after relabelling some of 
the summation indices) the discrete Euler-Poincare equations ([S]). □ 

In this context, a variation can be interpreted in two ways. In the case of unreduced 
fields, a variation of a field is a map 50 : TG such that b<^ij G T^,. .G. In the 
case of reduced fields, a variation of a field {wjj, f jj} is a pair of maps bu, bv -.V ^ TG 
such that bui^j G Tu^^G and G Ty-.G. In both cases, we demand that the variation 
is zero on the boundary of U: 



0, and buij = bvij = for all {i,j) G dU. (11) 



Here, the boundary dU is defined as the set of vertices such that {i,j) is a vertex 
of at least one face in U, and at least one face not in U (see [T6]). 



3.3 The reconstruction problem 

From theorem 13.41 we know that a solution (p : V ^ G of the discrete Euler-Lagrange 
equations gives rise to a reduced field (p : E ^ G, which has a natural interpretation 
as a fiat discrete connection in the sense of definition 12.51 
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To tackle the converse problem, we use the following consequence of proposition I2.6[ 
Recall that the holonomy mapping a; of a flat discrete connection u; is a map from V xV 
to G. 

Proposition 3.5. Let uj be a flat discrete G-connection with associated discrete holon- 
omy (2j : V X V ^ G. Then there exists a map cf) : V ^ G such that uj{xo,Xi) = 
(l){xo)~^(p{xi) . The map is unique up to left translation by an element of G. 

Proof. Choose an arbitrary vertex xq and a group element go, and define 0(xo) = go- 
Let Xi be any other vertex and put 0(xi) = goU}{xo,Xi). This map is well defined. □ 

Let if : E G he a. solution of the discrete Euler-Poincare equations ([8]). We now wish 
to reconstruct a solution (p of the original problem, such that = $ o 0. The map (p is 
provided by proposition 13.51 on the condition that is a flat connection. As soon as 
is not fiat, the holonomy ci) is path dependent, and no such (p can exist. Therefore, we 
have the following theorem. 

Theorem 3.6 (reconstruction). Let (p : E ^ G be a solution of the discrete Euler- 
Poincare equations Iji^. There exists a solution cp : V ^ G of the unreduced Euler- 
Lagrange equations ^ if and only if (p is flat. In that case, (p is uniquely determined 
up to left translation by an element of G. 

Remark 3.7. In some cases, the element go & G used in constructing the map (p is 
fixed by considering boundary or initial conditions on the solutions. o 

4 The Noether theorem 

In the case of continuous field theories, Noether 's theorem states that, for each contin- 
uous symmetry, there exists a conservation law (see [7] for an overview). Here, we will 
show that a similar theorem holds for discrete field theories with a continuous symme- 
try, and that the conservation law associated to the left G-invariance of the Lagrangian 
is equivalent to the Euler-Poincare equations. A similar theorem was proved in [5] in 
the case of Lagrangian reduction for continuous field theories. 

4.1 The Poincare-Cartan forms 

We begin by introducing a set of Poincare-Cartan forms that will be instrumental in 
formulating the Noether theorem. It should be remarked that the definition of the 
Poincare-Cartan forms here is a special case of a more general construction, elaborated 
in [20] and briefly outlined in section [61 

As in section [3l let L : = G x G x G — M be a G-invariant Lagrangian and consider 
the reduced Lagrangian / : G^ = G x G ^ R. Associated to L is a set of three one-forms. 
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called Poincare-Cartan forms, and defined as follows: 

with Vi G Tg^G, and similarly for 6*^^ and 9^^y Note that + 6*^^ + 6*^^ = dL. 

For the reduced Lagrangian /, the definition of the Poincare-Cartan forms is somewhat 
less direct. We put 

Ol):G^^g*(BQ*(BQ\ (12) 
{e[,){u,v),{^,,^2,^3)) = -{dl{;v),TR^{^,)) - {dl{u,-),TR,{^,)), 

where ^1,^2,^3 ^ 0- Alternatively, if t 1— h(t) is a curve in G such that h{0) = e and 
h{0) = 6'|-,^^ may be defined as 



(^(i)K^^),(ei,6,e3)> = ^i{h{t)-\h{t)-'v) 



i=0 



dt 

The remaining one-forms O''^^) ^(3) defined by 

(6,6,^3)) = (d/(-,t;),TL„(6)) (13) 

and 

(^(3)(«,^),(ei,6,e3)> = {dl{u,-),TL,{^,)). (14) 

The apparent asymmetry between (fT2|) and (fT3| [T^ is due to our definition of /. For 
a fully symmetric set of Poincare-Cartan forms, one should follow the prescriptions of 
EOl. 



The Lagrangian L is related to the reduced Lagrangian / as follows: L = / o ^. A 
straightforward computation (or the application of theorem 25 in [20]) shows us that 
the corresponding Poincare-Cartan forms are related in a similar way: 

Lemma 4.1. Let \E' be the bundle map {'^,'^) as defined above. Then 



Furthermore, a similar identity holds for the presymplectic forms fif := — d^J.-, and 



:= -d9^ ■ 



We recall that the puUback of 9[^^ by the bundle map ^ = (^, ^) : T^G f ^ (G x q) 
is defined as follows: 



((^*^j,))(^l, ^2, ^3), (^'1, V2, V3)) = (^j.)(^(^i, g2, ^3)), ^(^^1, V2, V3 
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Figure 3: Location of fi, f2, and fs. 
4.2 The unreduced Lagrangian 

The Lagrangian L : ^ R is assumed to be left G-invariant in the sense that 
L{ggi, 992, ggs) = L{gi, g2, gs) for all g in G. According to Noether's theorem, which 
we will prove in a moment, there is a conservation law associated to this symmetry. 

Let ^ be an element of g. Infinitesimal invariance of Lagrangian under the flow generated 
by ^ is expressed as 

{dL{g,,g2,gs),{^G{gi),^G{g2),^G{g3))) = 0, (15) 

where ^g, defined by C,G{g) = TRg{^), is the fundamental vector field associated to ^. 
Consider the functions : ^ M, z = 1, 2, 3, given by 

Jl{gi,g2,gz) = (^(^)(^i,^2,^3), (^g(^i),^g(^2),^g(^3))) ■ 

As J| is linear in ^, we can define a map : ^ g* by the prescription {J\C,) = J^, 
for all G g. An immediate consequence of the G-invariance is that 

J^ + J^ + J^ = 0. (16) 

Now, let there be given a map (f) : V —>■ G and consider the pull-back of each J* by 
ijj := (p X (p X (f). In this way, we obtain a map from F to g*, which can be identified, by 
means of the discrete Hodge star, with a map from V* to g*. In view of f fT6|) it turns 
out that we only need to consider two such maps, which we will denote by r]^^^ and 7]^y\ 
and which are given by 

r^W(r) = J\<P{X,),<P{X2),<P{X3)), and v'^y\r) = j3(0(xi), ^(xs), ^(xg)) 

for r G V*, and where f = (xi, X2, Xs, X4) is the face dual to r: *f = r. We further put 
^ <r/(^),^) and r]^^^ = {v^y\^). 

Before stating Noether's theorem, we recall the definition of the backward difference 
operators 6~ and S~: if ri, r2 and r^, with rj = *fj, are located as in figure [31 and 
/ : V"* — >■ M is a function, then 

(^.T/)(ri) = /(ri)-/(r2) and {6^ f){r,) = f{n) - f{r,). 
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These difference operators can be extended without difficulty to the case of vector- 
valued functions. 

Proposition 4.2 (Noether). Consider a G-invariant Lagrangian L : ^ M.. If cf) is 

a solution of the discrete field equations then the maps r]^^\r]^y^ : V* ^ Q* satisfy 
the following conservation law: 

s-v^""^ + s;v^'^ = 0. (17) 

Conversely, if (p is such that the associated mappings ri^^'> and r]^^^ satisfy [T7\ l, then (p 
is a solution of the discrete field equations. 

Proof. Consider points ri,r2,r3 of the dual mesh as in figure [3] and let be the dual 
face to rf. f j = *rj. Then we have that 

= Jimi)) - Jlim)) + Jimh)) - Jimfs)) 
= -j'M^i)) - J!m2)) - j'ims)), (18) 

where we have used (fT6ll . Recalling that ip is defined as ipif) = (0(^i)) 0(^2), 0(2^3)), 
we note that 

Therefore, we can rewrite the conservation law ( fTSi) as 

/Pit Pit f)T 

which vanishes if (p satisfies the field equations. Conversely, if the left hand side is zero, 
then is a solution of the field equations. □ 



4.3 The reduced Lagrangian 

Not only is Noether's theorem equivalent to the unreduced discrete field equations, 
it will turn out that it contains the Euler-Poincare equations as well. To show this, 
we start from the discrete conservation law as expressed by (fT8|) . We use the same 
notational conventions as in the preceding section and write 

^{^{U)) = {u,,v,), ^ = 1,2,3. 

Furthermore, we note that (V'(fi))i = (V'(f2))2 = ('?/'(f3))3; this unique element of G is 
denoted by g. By rewriting each of the three expressions in ( ITSl) in terms of the reduced 
Lagrangian / only, we obtain 

4mi)) = -{R:,dii;v,),v) - (i?;d/(«i,-),r7>, 
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where rj = Ad^-i^, as well as 

Jlm,)) = {Lldli-^^^)^v) and Jim,)) = (Ll^dliu,, ■ 
Putting all of these expressions together gives the following: 

6-vt\ri)+S;vfiri) = ( [R:,dl{;v,) - L:^dl{;v,)] + 

[R:^dliuu-)-L:^dl{us,-)],v). (19) 

As rj ranges over the whole of g, we conclude that the conservation law (ITTl) implies the 
discrete Euler-Poincare equation ([8]). 

5 Extending the Moser-Veselov approach 

In their seminal paper, Moser and Veselov [17j approached the problem of finding an 
integrable discretization of the rigid-body equations by embedding the group SO{n) 
into a linear space, namely gl{n). Somewhat later, Marsden, Pekarsky, and Shkoller 
[15] then developed a general procedure of Lagrangian reduction for discrete mechanical 
systems, and showed that the Moser-Veselov equations are equivalent to the discrete 
Lie-Poisson equations. 

Here, we intend to do the same thing for a fundamental model in field theory: that of 
harmonic mappings from into a Lie group G. We will show that it is possible to 
develop a Moser-Veselov type discretization of these field equations, provided that G 
is embedded in a hnear space. As could be expected, these discrete field equations are 
equivalent to the Euler-Poincare equations. 

In the continuous case, the harmonic mapping Lagrangian is given by 

c = \ (0" v., r v.) + 1 (r v„ r v.) , (20) 

where (■, ■) is the Killing form on g and where the subscript 'x' and 'y' denote partial 
differentiation with respect to that variable. For the sake of clarity, we will only treat 
the case of harmonic maps that take values in SO{n), embedded in gl{n), in which case 
the Killing form is just the trace. We stress that the entire theory can be generalized 
to the case of an arbitrary semi-simple group G, embedded in a linear space. 

Consider the rectangular lattice from section 12.11 and denote the lattice spacing by h. 
As usual, we denote the values of the field (j) on the lattice points by We discretize 
the reduced partial derivatives (p'^cpx and (p~^(py by writing them as follows: 

0~Vx ~ T(Pl+i,j{(Pi+id - <Pi,j) and (p''^(py ^ \4>lj+i{<Pi,j+i - (pij), 
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where (pjj is the transpose of the matrix (pij. Substituting this into (l20l) yields the 
following discrete Lagrangian (up to an unimportant constant): 

Cd = -^tr(0j.0i+i,j-) - ■^tr(0j'^.0ij+i). 

In order to ensure that (f>ij G SO{n), we need to impose the constraint that (pjj4>i,j = ^■ 
We are thus led to consider the following action: 

S{<P) = J2 (tr(C-0^+l,,) + tTi<|>y^,+l) - {^^A<ff,<l>^,J -I))), (21) 

where we have rescaled the Lagrange multipliers Aij to get rid of the factor — 
Note that Aij is a symmetric matrix. 

The field equations are obtained by demanding that S be stationary under arbitrary 
variations; they are given by 

(pi+ij + (pi-ij + (pij+i + (pij-i = (pijAij. (22) 

We multiply these equations by 4>Jj from the right, and use the fact that (f)ijAij(j)fj is 
a symmetric matrix in order to get rid of the Lagrange multipliers: 

By introducing the following quantities, 

m+ij = - and Uij+i = - ^ij^^j+i, 

the field equations can be rewritten as the following set of conservation laws: 

rrii+ij + Uij+i = niij + Uij. (23) 

Finally, let us introduce the discrete momenta Mij and Nij, defined as 

= (f)J_^ijmij(j)i-ij and Nij = (f)Jj_^nij(j)ij_i. 

The field equations governing the behaviour of these quantities are then easily deter- 
mined to be, on the one hand 



Mi J = ttij - ajj where a^j = Uj-ij, ^^4) 
= PiJ - where /3ij = vtj-i 



as well as, on the other hand, the counterpart of (123|) : 

M+i,, + N,j+^ = Ad^rMij + Ad^rNi^j. (25) 

The similarities with the Moser-Veselov equations for the discrete rigid body are obvious 
(compare with equation (4) in [H]). 
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Remark 5.1. It is now straightforward to see the equivalence between the Moser- 
Veselov and the Euler-Poincare equations. Indeed, starting from the reduced La- 
grangian /, put 

= Vij) and Nij+i = Rl^ Al{uij, •), 

which can be interpreted as a discrete Legendre transformation [20]. Furthermore, put 
aij = Ui-ij and Pij = fij-i. The Euler-Poincare equations (IHl) then reduce to the 
Moser-Veselov equations derived above. o 

5.1 Relation with the Euler-Poincare equations 

Instead of deriving the Euler-Poincare equations from the Moser-Veselov equations, we 
can also start directly from the discrete action sum fl?I]) and proceed as in section [31 

The unreduced field equations are given by fl22l) . We now multiply these equations from 
the left (rather than from the right as in the derivation of the Moser-Veselov equations) 
to obtain the following set of discrete Euler-Poincare equations: 

Uij + Vij + uf_^j + vfj_^ = Aij, (26) 

together with the integrability condition 

UijVi+ijUl-j^^v^-j = e. (27) 

Using then the symmetry of A, .,, we eliminate the multipliers Ajj to arrive at the 
following expression: 

If we view as a 0[(n)-valued discrete one-form in the sense of §212.21 . then the Euler- 
Poincare equations can be conveniently expressed using the discrete codifferential: 

5v = {^Vf. (28) 

or \5(p\ = 0, where [■] denotes the antisymmetric part of a matrix: {A} = |(A — A^). 
Many authors (see, for instance, [231 15]) express the continuum equations for harmonic 
mappings taking values in a Lie group G with bi-invariant metric as da = = 6a, 
where a is a g-valued one-form on the base space. Equations flTTj) and fl25|) are the 
discrete counterpart of these continuum equations. 

6 Generalizations 

6.1 Field theories on non-trivial base spaces 

Up until now, we considered discretizations of field theories with two independent vari- 
ables, i.e. where the base space is M^. From the beginning of the paper, we discretized 
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by considering the square mesh {V, E, F) in whose vertices have integer coordi- 
nates. Other discretisations of can be introduced as in |20| by noting that (V, E, F) 
is a planar graph; the theory of this paper can be readily generalized to the case of 
arbitrary planar graphs. 

Secondly, the generalization to an arbitrary base space X, not necessarily M^, is straight- 
forward under the assumption that a simplicial complex can be embedded in X. Let 
us now briefly discuss this case. 

As before, unreduced and reduced fields are maps associating a group element to each 
0-dimensional and 1-dimensional simplex, respectively. Other generalizations are also 
possible, for instance where one considers fields that associate group elements to higher- 
dimensional simplices, but these types of fields do not arise in Euler-Poincare reduction 
and cannot be treated with the Lie groupoid framework of section 16.21 For an intro- 
duction to such field theories, see [22] . 

The theory of discrete differential forms starting from a simplicial complex and its dual 
is developed to great detail in [8], [lO] for spaces of arbitrary dimension and topology. 
To study discrete Lagrangian field theories, one now has to focus on the Cartesian 
products and G^", where n -|- 1 is the dimension of the base space. The usual 

procedures of discrete Lagrangian field theories and discrete reduction of the previous 
sections carry through to this case without significant changes. 

A significant change occurs when the base space is not simply connected. In that case, 
proposition 12.61 is no longer valid. Consequently, a reduced field then no longer induces 
a uniquely determined unreduced field. This is similar to the continuous case; see [Sj. 

6.2 Discrete field theories on Lie groupoids 

For the better part of this paper, we have considered discrete field theories associating 
group elements to either the vertices (unreduced fields), or to the edges of the mesh 
(reduced fields). It turns out that both are special cases of discrete field theories taking 
values in Lie groupoids. This point of view was introduced in [20], to which we refer for 
a detailed overview. In this section, we will briefly show how Lie groupoid field theories 
provide a unified framework for Euler-Poincare reduction. We will focus mostly on a 
Lie groupoid version of theorem 13.41 but note that other constructions, such as the 
Noether theorem and the Poincare forms of section HI also have their counterpart in 
the Lie groupoid framework. 

For more information on Lie groupoids, as well as their role in discrete mechanics, the 
reader is referred to [I3l EH HI] . 

A Lie groupoid is a set G with a partial multiplication m, a subset Q of G whose 
elements are called identities, two submersions a,(3 : G ^ Q (called source and target 
maps respectively), which both equal the identity on Q, and an inversion mapping 
i : G ^ G. A pair of Lie groupoid elements {g, h) is said to be composahle if the 
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multiplication m{g,h) is defined; the set of composable pairs will be denoted by G2. 
We will denote the multiplication m{g,h) by gh and the inversion i{g) by g^^. In 
addition, these data must satisfy the following properties, for all g,h,k E G: 

1. the pair {g, h) is composable if and only if [3{g) = a{h), and then a{gh) = a{g) 
and p{gh) = (3{h)- 

2. if either {gh)k or g{hk) exists, then both do, and they are equal; 

3. a{g) and /3{g) satisfy a{g)g = g and gl3{g) = g; 

4. the inversion satisfies g~^g = (3{g) and gg~^ = (y{g). 

A Lie group G can be considered as a Lie groupoid over a singleton {e}: the source and 
target maps a and (3 map any group element onto e and the multiplication is defined 
everywhere. 

Another example is the pair groupoid Q xQ, where Q is a manifold. The pair groupoid 
is a Lie groupoid over Q, and the source and target mappings are defined as follows: 
(^{Qo^Qi) = Qo: and /3(go,5'i) = Qi- The multiplication of composable elements is then 
given by (go, ^2) = (go, q2)- The set of units in Q x Q is {{q, q) : q E Q}. 

A groupoid morphism is a pair of maps (f) : G ^ G' and f '■ Q ^ Q' satisfying 
a' o (f) = f o a, 13' o (j) = f o p and such that (f){gh) = (j){g)(j){h) whenever {g, h) is 
composable. Note that {4>{g), (t>{h)) is a composable pair whenever ((?, h) is composable. 

6.2.1 Discrete fields 

The set of edges is a subset of the pair groupoid V x V , but it is not a groupoid 
in itself because the multiplication of two elements of E is not necessarily again an 
element of E. Nevertheless, it turns out that much can be gained from considering E 
as a "local groupoid" in the sense of [201, rather than just as a subset of x V^. By 
a "local groupoid", we mean here that, even if the multiplication is not defined in E, 
one can still define composable edges as pairs of edges ei,e2, such that (3{z2) = «(ei), 
where a{xo,xi) = xq and P{xq,xi) = xi, i.e. a and (3 are the usual source and target 
mappings of the pair groupoid V x V, restricted to E. Hence, E is a. groupoid in all 
aspects but one, the multiplication. 

With the conventions introduced above, we have for e = {xo,xi) that = {xi,xq). 
In addition to the mesh (V, E) in M^, we now consider an arbitrary Lie groupoid F 
over a manifold Q. The idea is to define discrete fields as mappings from the "local 
groupoid" E to the Lie groupoid F. In particular, if (y9 is a discrete field, then ip maps 
composable edges in E to composable elements in F, and, secondly, if c is any edge in 
E, then ip{t~^) = This is worked out in more detail in the following definition: 

Definition 6.1. discrete field is a pair if = {<f{o), where ifi^o) is a map from V 

to Q and ^9(1) is a map from E to T such that 
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1. a(v3(i)(x,?/)) = v2(o)(a;) and f3{<^(^i){x,y)) = ip(^o){y); 

2. for each {x,y) G E, V5(i)(z/,a;) = [(p(^i){x,y)]'^ ; 

3. for each x eV, {p(i){x,x) = V5{o)(a^) G Q- 

Even though we are working with an object E which is not quite a groupoid, it turns 
out that any discrete field ip : E ^ T may be extended to an actual groupoid morphism 
from V xV to r. This can be shown quite easily; for a proof, we refer to proposition 7 
in [20] . Note that proposition 12.61 is a special case of this extension property. 

Example 6.2. Consider a discrete field (^9(0), taking values in the pair groupoid 
Q X Q. Because of the conditions in definition 16. the map ip(^i) : V x V ^ Q x Q, 
is equal to (y9(o) x (/}(o)- In this case, the discrete field is completely specified once we 
are given (y9(o) : V ^ Q. For Q = G, this corresponds to the class of unreduced fields 
studied in section [31 o 

Example 6.3. Consider now the case where the Lie groupoid F is a Lie group G and 
let ((y9(o), V2(i)) be a discrete field. As ip(^o){x) = e for all x E V, a. discrete field can 
be identified with a map : V x V ^ G. As in section 13.31 this is a fiat discrete 
connection, or, equivalently, a reduced field. o 

6.2.2 Reduction 

The previous two examples show that the Lie groupoid framework encompasses both 
the case of unreduced and reduced fields (when F = G'xG'orF = G', respectively). 
The following theorem shows that the Euler-Poincare reduction of one Lie groupoid 
field theory yields another one: 

Theorem 6.4 (see [20]). Let T' be a Lie groupoid over a manifold Q' and consider a 
morphism ($, /) : (F, Q) — > (F', Q'). Furthermore, let L' : G'^ — >■ R 6e a Lagrangian on 
G'^ and consider the induced Lagrangian L = L' on G^ , where \& : G'^ ^ G'^ is the 
map associated to $. 

A morphism : x V" ^ F will satisfy the discrete field equations for L if the induced 
morphism ^ o cj) : V x V V satisfies the discrete field equations for L' . 

Lemma 6.5. // the morphism ($, /) in theorem 6.4 is a submersion, then a discrete 



field (p : V ^ V T will satisfy the discrete field equations for L if and only if the 
induced field ^ o (p : V x V ^ T' satisfies the field equations for L' . 

Proof: The proof follows that of Corollary 4.7 in [T^. o 

A few remarks are in order here. The manifold G^ consists of sequences of k composable 
elements {gi,g2, ■ ■ ■ , Qk) in F, such that (71 ■ (72 ■ ■ ■ fi'fc is a unit: 

G'' = {(^1, ^2, • • • , gk) e : a{gi+i) = (3{gi), a{gi) = (3{gk) and 

gi- g2---gk = ea{g^)}- 



21 



The manifold is the manifold where the Lagrangian is defined, and as such it is the 
generalization to the case of Lie groupoids of the Cartesian products G x G x G and 
GxG used in section I3TT1 It is straightforward to check that, for k = 3,G^ = GxGxG 
when r = G X G, and = G X G for r = G. 

The map $ : G x G ^ G, introduced in section [STTl and defined as ^{g,g') = g~^g', 
is a Lie groupoid morphism. When we put F = G x G and V = G, and use $ in 
theorem I6.4[ we obtain the reduction theorem I3.4[ In contrast to discrete mechanics, 
it is not very likely that a reconstruction theorem can be proved in the context of Lie 
groupoids; in the case of Lie groups, for instance, there is the obstruction of flatness. 

7 Concluding remarks 

In this paper, we have shown that the well-known concepts of symmetry and Euler- 
Poincare reduction can be extended by very modest means to the case of discrete field 
theories. It turned out that many of the constructions known from continuum field 
theories have a natural discrete counterpart. However, much remains to be done. 

Throughout the text, we remarked that almost all of our constructions can be extended 
to the case of field theories taking values in an arbitrary Lie groupoid F. Such a 
framework might be useful in the case of general Lagrangian reduction. In addition, 
as we pointed out before, many of the definitions introduced here gain much in clarity 
when rephrased in this more general language. 

As a second point of interest, we are also interested in the concept of multisymplectic- 
ity. Because of their variational nature, it is to be expected that the Euler-Poincare 
equations are multisymplectic in the sense of H]. However, as the focus of this 
article was on theoretical developments rather than on the construction of practical 
integration schemes, we have chosen to leave this topic for future work. 
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